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Abstract 

The roles of natural convection In the melt and the shape of the melt/solld 
Interface on radial dopant segregation are analyzed for a prototype of vertical 
Bridgman crystal growth system by finite element methods that solve simultaneously 
for the velocity field In the melt, the shape of the solidification Isotherm, 
and the temperature distribution In both phases. Results are presented for 
crystal and melt with thermophysical properties similar to those of gallium- 
doped germanium In Bridgman configurations with melt below (thermally destabi- 
lizing] and above (stabilizing) the crystal. Steady axlsymmetrlc flows are 
classified according to Rayleigh number as either being nearly the growth 
velocity, having a weak cellular structure or having large amplitude cellular 
convection. The flows In the two Bridgman configurations are driven by different 
temperature gradients and are In opposite directions. Finite element calculations 
for the transport of a dilute dopant by these flow fields reveal radial segrega- 
tion levels as large as sixty percent of the mean concentration. Segregation 
Is found most severe at an Intermediate value of Raylelnh number above which 
the dopant distribution along the Interface levels as the Intensity of the 
flow Increases. The complexity of the concentration field coupled with 
calculations of effective segregation coefficient show the coarseness of 
the usual diffusion-layer approximation for describing the dopant distribution 
adjacent to the crystal. The length of the gradient zone and the ratio of 
thermal conductivities between melt and crystal are Identified as additional 
critical parameters for setting the degree of radial segregation. The predic- 
tions of melt/solid Interface shape and axial and radial dopant segregation 
will form the basis for the first comparison between detailed fluid mechanical 
models of melt crystal growth and actual experiments. . 
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1. INTRODUCTION 

The compositional uniformity of dopants and Impurities In semi- 
conductor crystals grown from the melt depends strongly on the pattern 
and Intensity of flow In the melt and on the shape of the melt/sol Id 
Interface. In processes such as the Czochralskl and floating zone 
methods for producing silicon and the vertical Bridgman configuration 
studied here, the unequal partitioning of dopant between melt and crystal 
during solidification causes a concentration gradient perpendicular to the 
melt/solld Interface. When the melt/solld Interface Is planar and the melt 
Is quiescent except for motion caused by crystal growth the concentration 
field decays exponentially with distance Into the melt. Curvature of the 
solidification Interface and convection In the melt both change the concentra- 
tion distribution along the Interface and alter the dopant level In the 
crystal . 

The effect of convection on segregation In the growth direction was 
pointed out by the Burton et al. analysis [1-2] of the effect of crystal 
rotation on mass transfer. When the crystal diameter Is large and the 
velocity field Is similar along the radius of the melt/solld Interface 
convection only alters the concentration field perpendicular to the Inter- 
face. Burtln et al. analyzed this case and developed an expression for an 
effective segregation coefficient and Introduced the notion of an axial 
"boundary layer" thickness for diffusion controlled mass transfer. Others 
[4-8] have applied the Idea of a uniform diffusion-layer or stagnant-film 
that Is masking a growing crystal from a well mixed bulk and have determined 


so-called boundary-layer (more appropriately called diffusion-layer) thick- 
nesses without any picture for the fluid motion In the melt. Other than 
Burtln, Prim and Sllchter's original analysis of rotating flows, few papers 
dealing with melt crystal growth [9-10] have addressed the exact coupling 
between the fluid flow and mass transfer In a consistent way. 

In many growth configurations, the concept of a uniform diffusion -layer 
yields an over-simplified view of the role of convection In dopant segregation 
[11]. When the flow Is laminar and cellular, as Is the case for many 
small-scale growth systems and for reduced gravity experiments, convective 
mass transfer Is uneven along the surface of the crystal and significant 
radial segregation results. The coupling between moderate convection and 
radial dopant segregation has been thoroughly analyzed for the rotational ly- 
drlven flows In small-scale floating zones [10, 12, 13] and has been demonstrated 
for buoyancy-driven convection under microgravity conditions [14], All 
these studies were of model crystal growth systems with planar solidifica- 
tion Interfaces. 

Corlell and Sekerka ([15]; also see ref. 16) have demonstrated that sig- 
nificant radial segregation occurs In systems without convection tangent to the 
crystal surface when the radius of curvature of this Interface Is the same 
or less than the length scale of the concentration gradient adjacent to the 
Interface. Curvature -Induced segregation has been shown to be an Important 
contribution to dopant Inhomogeneities In capillary growth systems [17-19] 
and may be Important In other small-scale crystal growth experiments In 
which convection In the melt has been suppressed. 

The purpose of this paper Is to present theoretical results from 


- 3 - 


computer-aided analysis of the Interactions of natural convection and melt/ 
solid Interface shape In setting dopant distributions In crystals grown In 
a vertical Bridgman system. The prototype growth system considered here 
Is shown In Figure 1 and consists of crystal and melt contained In a cylin- 
drical ampoule of radius R and length L pulled slowly through a vertically 
aligned furnace with hot and cold Isothermal zones. The Isothermal zones 
are separated by an adiabatic region designed to promote steep axial tempera- 
ture gradients and to maintain a flat solidification Interface. Heat transfer 
In the Bridgman system without convection Is well understood. A host of 
one- and two-dimensional heat transfer models [20-27] have been used to analyze 
the magnitudes of axial and radial temperature gradients In the gradient zone 
and the shape of the melt/sol Id Interface. No previous studies have con- 
sidered the effect of convection on melt/sol Id Interface shape. 

Convection In the Bridgman system Is driven by buoyancy differences 
Induced by radial and axial temperature gradients and may be either laminar, 
perfectly time-periodic or turbulent, depending on the magnitude of the 
temperature gradient, the size of the ampoule and Its orientation with 
res. ct to the direction of gravity. When the melt. Is positioned below the 
crystal the axial temeperature gradient places the hottest (and less dense) 
melt on the bottom of the ampoule and Is destabilizing. Kim et al. [5] 
used dopant strletlons In gallium-doped germanium crystals to Identify 
melt lengths at which each of the three forms of convection were present. 

Others [8, 28] have oriented the melt above the crystal In order to produce 
a stably-stratified axial density gradient. Convection still exists In this 
configuration; radial temperature gradients Introduced by changes In lateral 
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heat transfer conditions cause lateral density variations that can drive 
laminar and even turbulent convection, as demonstrated In [29]. 

The adverse effects of large amplitude convection, especially time- 

dependent flows, on dopant segregation In the crystal are so severe 

that solidification experiments have been performed In outer space [30, 31] 

to test the feasibility of crystal growth In a reduced gravity environment. 

Even these experiments are not convectlonless because time-averaged accelera- 

-4 

tlons with magnitude as large as 10 the gravitational acceleration on earth 
still exist on board most spacecraft. As discussed below, the differences 
In space experiments may cause larger concentration variations than seen 
In crystals grown on earth.- 

The calculations described here are of the axlsymmetrlc steady state 
velocity, temperature, and concentration fields and melt/sol Id Interface 
shapes for the model of the vertical Bridgman system described In Section 2. 
The mathematical free-boundary problem that describes these variables 
Is solved by a newly developed [31, 32] finite-element/Newton technique 
that computes simultaneously the shape of the solidification Interface, 
the velocity and pressure fields In the melt and the temperature distri- 
bution In both melt and crystal. The increased efficiency of this algorithm 
over any other developed before [31] makes possible the wide range of 
calculations presented here. Also, as described In [32], the use of 
Newton's method for solution of the full set of nonlinear algebraic equa- 
tions makes available powerful algorithms for detecting either 
the loss of solution or the existence of multiple solutions 
to the free-boundary problem represented by the finite-element equations. 


Flow fields calculated with the finite-element/Newton algorithm a^e used 
In a separate calculation of the distribution of a dilute dopant within 
the melt and crystal. All the numerical methods are presented In detail 
elsewhere and so are only sketched In Section 3. 

Flow fields and concentration distributions In the vertically stabilized 
(melt on top) configuration are presented In Section 4 along witr; sensitivity 
analysis for the length of the melt, the size of the gradient zone and the 
ratio of thermal conductivities between melt and solid. The usefulness of 
the concept of a diffusion-layer thickness for correlating the segregation 
results Is addressed In Section 4.2. The flow patterns and segregation 
results predicted for this growth configuration are compared In Section 5 
to results for the vertically destabilized (melt on bottom) arrangement. 

The bulk of the calculations reported here are for a melt and crystal with 
thermophysical properties similar to those of the gallium-doped germanium 
system used In several small-scale experiments [5, 8, 29]. The predictions 
of melt/solid Interface shape, radial dopant profiles and effective axial 
segregation coefficients are used to establish the link between the 
calculations and experimentally observable quantities. 
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2. PROTOTYPE OF VERTICAL BRIDGMAN GROWTH SYSTEM 

The finite length of a real ampoule In a Bridgman furnace causes the 
heat and mass transfer Inside the sample to be time dependent for any 
nonzero translation rate. These transients are caused by the steady 
decrease In the ratio of melt to crystal In the ampoule and are slow for 
the growth rates typical for semiconductor materials. We replace the true 
unsteady problem with a prototype steady state process viewed from a 
stationary reference frame and described In the cylindrical polar coordinate 
system shown In Figure 1. The translation of the ampoule Is accounted for 
by supplying melt Into the ampoule at z ■ 0 with velocity v 2 * and remov- 
ing crystal at the other end of the ampoule at the rate V $ = (Pj/p^V* * cjV £ 

For long ampoules and melts with low Prandtl numbers, the transient 
effects on heat transfer are small [26] and the changes In the fluid 
mechanics are accurately accounted for by the steady state model with 
calculations for a succession of ampoule positions Inside the furnace, 
each of which corresponds to. a different time In the actual growth experi- 
ment. Changes In the length of the melt caused by translating the ampoule 
lead to transients In dopant transfer that force the average concentration 
of dopant over the cross section of the crystal to change with time. Our 
hypothetical steady state model for mass transfer Ignores this effect. 

The prediction of effective segregation coefficients caused by convection 
from the results of the steady state model Is addressed In Section 4.2. 

The field equations governing heat transfer In both phases and fluid 
flow In t'le melt are put In dimensionless form by scaling lengths with the 
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height of the ampoule L, velocity v(r,z) with o^/L where Is the thermal 

dlffuslvlty In the liquid (see the list of thermophyslcal properties In 

Table I), pressure p(r,z) with p A a £ /L and temperature e(r.z) s (T(r,z) - T c )/ 

(T h - T c ). In this last definition T(r,z) Is the dimensional temperature 

field and T„ and T u are the temperatures of the hot and cold regions of the 
c n 

furnace. In units of ampoule length, the dimensions of the ampoule and 
furnace are the aspect ratio of the ampoule A = R/L, the length of the hot 
zone L h 2 L^/L and the length of the gradient region L g = L g /L. These dimen- 
sionless scalings, especially the choice of L as the length scale, have been 
made for convenience and may not give dimensionless groups that are the most 
appropriate for the physics controlling the buoyancy-driven flow and associated 
dopant segregation. The formulation more meaningful scaling Is addressed 
In Section 4.5. 

The location of the melt/solld Interface Is represented as z « h(r), 

0 < r < A , and the unit vectors everywhere normal N and tanoent t to the 

“ mm. 

Interface are 


*i * h r e r 

Ns— r * r 


(T“h7) 1/J • 
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5r 4 h r!i 


(1 ♦ h r ) 




0 ) 


where h r 2 dh/dr and the set (e p , e 0 , e,l are unit vectors for the cylindrical 
polar coordinate system. 

Field variables In the melt (0 <^ < h(r), 0 <. r < A) are governed by 
the axlsymnetrlc and dimensionless form of the Bousslnesq equations [32]: 


7 • v ■ 0 


( 2 ) 
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V • Vv • -Vp ♦ p r7 2 v + RaPr6e z . ^ 

v • 70 ■ 7*0 , (4) 

m 

where V = e p 3/3r ♦ 3/3 r Is the gradient operator In cylindrical coor- 

dinates. The Prandtl Pr and Rayleigh numbers Ra appearing In eq. (3) are 
defined as Pr = v/o ^ and Ra £ fig(T h -T c )L 3 /a,jV where the symbols are defined 
In Table 1L The velocity field In the crystal Is uniform at speed V $ and 
the temperature field In the crystal \[r) £ z ^ 1, O^r^A) Is governed 
by 

Pe e, • ?e - yV 2 e , (5) 

where y £ <* 5 /ct f Is the ratio of thermal dlffuslvltles In solid and melt and 
Pe £ VjL/cij, Is the Peclet number for convective heat transfer. 

The shape of the melt/solld Interface h(r) Is set by the conditions 
for the equilibrium meltlnq temperature T„ and the balance of conductive 
heat fluxes across the Interface with the latent heat released there. In 
dimensionless form, these Interface conditions are 

6(r,h(r)) 3 (T m - T c )/(T h - T £ ) , (6) 

(N • ve) t • K (N • 78). ■ SP* (N*e ) , (7) 

where S £ AH^/o^Cp " T c^ is the Stefan number and subscripts In 
JL 

eq. (7) denote whether the Interfaclal flux Is evaluated In the crystal 
(s) or the melt (i). 


Boundary conditions on velocity at the Interface Insure no slip tangen- 
tial to the crystal and the solidification of melt at a rate proportional to 
the growth rate V 8 and the difference between the densities a s (Pj/o^) 
of the two materials: 

v • t ■ Pe(e, • t) , o(vN) ■ Pe(e «N) . (8) 

« « *4 « m — *4 

The top surface and sidewall of the ampoule are assumed to be no slip 
surfaces so that v p ■ 0 and v 2 ■ aPe In the melt. The symmetry boundary 
conditions at the axis of the cylinder are 

v_ • 3v ,/3r ■ 36/3r *0, r ■ 0, 0 < 2 < 1 . (9) 

r z “ — 

The thermal boundary conditions along the ampoule are modelled after 
the 1deall2ed situation considered by Fu and Wilcox [2] In which the 
ampoule has negligible thermal mass, the adiabatic region Is a perfect 
Insulator and the heat transfer rates between the hot and cold regions 
of the furnace are so large that the ampoule has the temperature of the 
surrounding furnace. In addition, the ends of the ampoule are assumed 
to be perfectly Insulated. None of these assumptions Is totally realistic 
for small scale experimental systems. Actual ampoules provide resistance 
to radial heat transfer and conduct heat axially; this latter effect modifies 
axial temperature gradients [26] and causes radial temperature gradients 
near the melt/solld Interface [27]. The mathematical forms of the thermal 
boundary conditions along the ampoule wall are stated In [32] and are not 
repeated here. 
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The differential mass balance equation appropriate for dilute dopants 
Is 

(Sc/Pr) V • (cv) ■ V*c , (10) 

m 

where Sc s v/P Is the Schmidt number. The boundary conditions used for 
solving eq. (10) are 

N • 7c ■ (PeSc/Pr)(e_ • N) (1 - k)c, 0 < r < A, 2 ■ h(r) , (11) 

3c/3z ■ (PeSc/Pr) (c-1 ) , 0<^r^A , z-0 , (12) 

3c/3r * 0 , r**0»A , O^zO. (13) 

Equations (11) and (12) express conservation of mass at the melt/sol Id Inter- 
face and the fictitious "Inlet" at the melt end (z*0) of the ampoule, res- 
pectively. Equation (13) Is the no flux condltlorf valid at the centerline 
and sidewall of the ampoule. Once the velocity field In the melt and the 
shape of the melt/solld Interface are determined from the two-phase natural 
convection problem described above, equations (10-13) reduce to a linear 
set to be solved for the concentration distribution throughout the melt. 
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3. FINITE ELEMENT ANALYSIS 

The finite element/Newton technique developed In references [31-32] for 
solving steady solldlflcetlon problems has been applied to the model of the 
vertical Bridgman system described by eqs. (1-9). First, finite element 
approximations to the velocities, pressure, temperature and melt/sol Id 
Interface shape are combined with Galcrkln's method to reduce the original 
partial differential equations and boundary conditions to a set of nonlinear 
algebraic equations. For an approximate location of the solidification Iso- 
therm h(r), melt and solid are divided Into N p radial and N z axial elements; 
half of the axial elements are placed In each phase. A sample 4 x 16 (N f ■ 4, 

N z • 16) mesh Is shown as Figure 2. The field variables are approximated 
by expansions of standard finite element basis functions and coefficients 
which are to be determined. We use mixed Interpolation [33] for the Bousslnesq 
equations with continuous bilinear polynomials for Interpolating pressure and 
biquadratic ones for approximating temperature and velocities. The shape of 
the melt/solld Interface h(r) Is expanded In terms of Hermlte cubic functions. 

The f1?ld equations (2, 3, 5) are reduced by Galcrkln's method to algebraic 
residual equations In terms of the unknown coefficients for each variable. 

The energy flux condition eq. (7) Is Incorporated Into the residual equations 
for the energy equation (5) and the condition for the melting point Isotherm 
eq. (6) along the solidification Interface Is distinguished as the equation 
for determining the Interface location. A separate set of Interface residual 
equations are formed from eq. (6) by applying Galcrkln's method using the 
Hermlte cubic polynomials as weighting functions; see [31]. The entire set 
of nonlinear algebraic equations for the field variables and Interface 
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ship* are solved by Newton Iteration as described In [30, 31]. Computer- 
aided methods for detecting and tracking multiple solutions to the algebraic 
equations are also described In [32]. 

The equations for ness conservation of a dopant (10-13) are also 
solved by the Galerkln finite element method with the concentration field 
c(r,z) approximated by an expansion of unknown coefficients {c^} and blouadratlc 
basis functions {t 1 } 

c(r,z) ■£ c 1 ♦ 1 (r.t) , (14) 

1-1 

where Is the total number of basis functions In the expansion. Linear 
algebraic equations result from applylnq Galerkln* s method to eq. (10) and 
using the divergence theorem to Incorporate the boundary conditions 
01-13): N 

E Cj ■= b, . (15) 

1-1 


where 


A 1 j 3 f [♦’ V • (Vt 1 • V» J )] iV 

\ 

*/„ P.(1-k)(. -H)*VdA ♦ / P*«VdA , 

V I * Z * • / P 0 

t| * J ♦’ P. dA 
C 0 


06 ) 


(17) 


The nomenclature for the boundaries Interface Vj and Inlet and the melt 
volume V f Is shown on Figure 2. The equation set (15) Is solved by Gaussian 
elimination and the resulting concentration fields art plotted as contours 
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of c(r,z). 

Although eqs. (10-13) yield a linear problem, the solution Is by no means 
simple. The high velocities caused by natural convection and the low diffusl- 
vlties of dopants In melts lead to convection dominated dopant transfer. This 
Is clearly seen by computing a Peclet number for mass transfer In the melt 
as Pe n = V*L/tf, where V* Is the velocity scale for the convection. For 
natural convection driven by the vertical temperature gradient in the 
destabilized configuration, v* Is Ra^a^/L and Pe m = Ra^ 2 Sc/Pr * 1 x 10 5 

4 

for a Rayleigh number of 1 x 1C . Solution of the convection-diffusion 
problem (10-13) at this level of convection Is one of the outstanding 
problems In numberlcal analysis [34, 35], We have adopted two types 
of adaptive spacing of the finite element mesh to account for the steep 
concentration gradients that result from large Pe m . First, each element 
of the finite element mesh for velocity and temperature Is sub-divided 
Into smaller sections for approximating c(r,z). Secondly, these small 
elements are not equally spaced but are graded to be smallest near the 
boundaries of the melt where the concentration varies most rapidly. 


The stream function ip(r ,z) for each velocity field Is computed by 
solving the linear equation 


1 3 2 * 1 to 1 3 . 3v r 3v z 

7 77 *7 * 7 1 ? * ~ 


(18) 


by the Galerkln finite-element method with ^- ( r , z ) represented In a biquadratic 
basis, just as for the concentration field eq. (14). We display flow fields 
In terms of streamlines plotted as contours of ij/(r,z). 
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4. VERTICALLY STABILIZED CONFIGURATION 

Th« large number of dimensionless thermophyslcal and geometrical para- 
meters that must be set for any one simulation of the Bridgman system and 
the cost of each computation have forced us to concentrate on parameters 
similar to those for the gallium-doped germanium experiments underway In 
the Material Science Department at MIT [36]. The appropriate values of 
thermophyslcal parameters are listed In Table II and the corresponding 
dimensionless groups are given In Table I. The radius of the ampoule 
has been taken as 0.5 cm and the growth rate V $ has been set at 16 ym/sec. 

We first consider convection and segregation In the stabilized configuration 
(melt above w.'ystal) with the ampoule four times longer than Its radius, 
l.e. , A * 0.25. 

This configuration and the values of the dimensionless groups listed 
In Table II are the basis of the calculations presented In Sections 4.1 - 4.3. 
The sensitivity of the results to the ratio of thermal conductivities and 
to the design of ampoule and furnace are considered In Sections 4.4 - 4.5. 

4.1 Temperature and Velocity Fields 

Calculations of large amplitude natural convection In melts with low 
Prandtl number are extremely difficult and the accuracy of the numerical approdma- 
tlons can only be systematically checked by refinement of the finite element 
mesh. We have used grids of 4 x 16, 6 x 24, and 8 x 32 elements to compute 
the flows In the base case described above for Rayleigh numbers ranging from 
10 to 5 x 10 6 . Streamlines and Isotherms computed with the 8 x 32 mesh 
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are shown on Figures 3 and 5 for the test case. With this mesh, 239C nonlinear 
equations were solved by Newton's method; and each Newton Iteration required 
20 minutes of computer time In the Honeywell 6180 at MIT. 

Three distinct types of flow patterns were observed. At low Rayleigh 
numbers (Ra < 10^) the streamlines were rectilinear and only slightly distorted 
by buoyancy forces. For Intermediate values of Ra (10^ <_ Ra £ 10*) a 
cellular flow developed, which was driven by the radial temperature gradients 
established by the mismatch In thermal boundary conditions between the adiabatic 
and hot zones. The flow moved upward along the sidewall and downward at the 
centerlines of the ampoule. The center of the cell was located slightly 
above the gradient zone and migrated downward and toward the sidewall with 
Increasing Ra. Increasing the Rayleigh number to 2.6 x 10® led to the devel- 
opment of a weak secondary cell adjacent to the melt/sol Id Interface, as 
shown In Figures 3a and 3h; the onset of the multi-cellular flows marks 
the start of the third type of flow pattern. The motion In the secondary 
cell next to the Interface Is In the opposite direction to the main cell 
and leads to qualitatively different radial segregation than the flows 
with a single cell which exist for Ra less than 1 x 10®. This point 
Is discussed further below. 

The three types of flows found with varying Rayleigh number are 
also distinguishable on the plot of the circular of the primary cell 
shown as Figure 4. One Interesting feature of tills plot Is that a part of 
the Intermediate region can be Identified where \p |nax Is linear In Rayleigh 
number. This suggests that an exact linear analysis, such as performed 
for rotating flows by Harriott and Brown [12] and for buoyancy-driven 


- 16 - 


flows In simpler geometries by Dressier [37], Is applicable. Unfortunately, the 
temperature field In the melt, even In the absence of convection, Is suf- 
ficiently complicated that closed-form calculations are Impractical. 

The temperature fields for several values of Rayleigh number are displayed 
In Figure 5 and show the effect of the low Prandtl number in the pallium-doped 
germanium system. For Rayleigh numbers between zero and 1 x 10 , the thermal 
fields in both the melt and crystal are essentially the same as the field 
calculated without convection (Ra « 0) and are similar to results of Fu and 
Wilcox [21]. Increasing Ra above 1 x 10 4 caused the Isotherms along the 
axis of the melt to compress toward the melt/solid Interface by the downward 
fluid motion; by Ra * 5 x 10® the shape of several Isotherms farthest from 
the Interface has Inverted from convex to concave at the center of the melt. 

The Isotherms In the crystal were unchanged by changing Ra. Also, the 
large portion of the crystal at the uniform temperature 6 * 1 pointed out 
that the length of ampoule In the cold portion of the furnace was sufficient 
to guarantee that the position of the end of the ampoule was unimportant. 

The shape of the melt/solid Interface Is given by the Isotherm e * 0.5 
In Figure 5 and Is plotted separately on Figure 6. As Indicated by the 
temperature field the shape of the "Interface Is unchanged by convection 
for Ra between 0 and 10^. For higher Rayleigh numbers, the hot melt moving 
down the axis of the ampoule drove the melt/solid Interface deeper Into the 
adiabatic region. The changes In Interface shape caused by convection are 
not large; even for Ra ■ 5 x 10® the deflection of the Interface Is only 
six percent of Its mean location. 
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The fine finite element mesh (8 x 32 elements) was necessary for 
quantitatively accurate calculations over the ranqe 0 ^ Ra ^ 1 x 10® and 
for qualitatively accurate results above Ra * 1 x 10®. The accuracy of the 
calculations with the three different finite element grids was assessed 
by comparing results for and by computing the degree of closure of 
the overall heat balance as a function of Rayleigh number. These results 
are tabulated In Table III. The heat balance was computed as 

AQ s f (n*ve)dA/Nuj , (19) 

where the Integral was computed over all boundaries V M surrounding the 
melt and fhij was the Nusselt number evaluted at the solidification Interface 

(N • V6)dA . (20) 

I 

Calculations with all three meshes gave heat balances that closed to within 

one percent of the Interface flux and 4 ^ values that were within one percent 

4 

of the value for the 8 x 32 mesh when Ra was less than 1 x 10 . The errors 
In AQ and if/^ for the coarser meshes were larger at larger values of Ra. 

A more serious problem with the coarse mesh occurred for Rayleigh numbers 
greater than 1 x 10®. For both the 4 x 16 and 6 x 24 grids, sequences of 
calculations In Increasing Ra failed to converge for values of Ra between 
1 x 10® and 3 x 10®. Techniques for detecting multiple steady solutions of 
algebraic equations described In another publication [32] were used tt 
trace the families of steady solution through turns In Ra. The structure of 
these solutions Is represented on Fig. 7 by plotting the Nusselt number Nuj 
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as a function of Ra. For the 4 x 16 and 6 x 24 grids, three steady solutions 
to the algebraic solutions were found for several ranges of Rayleigh number. 

These multiple solutions are bogus; refining the grid to 8 x 32 resulted 
In the single flow family described above and also represented on Fig. 7. The 
reason for the failure of the coarse mesh calculations Is seen from the form 
of velocity field for Ra * 2.6 x 10 6 . The small secondary cells that are 
created alone the ampoule wall and at the melt/solid Interface were not 
resolved by the 4 x 16 and 6 x 24 grids. Clearly, coarse calculations of 
large amplitude convection for low Prandtl number melts can give qualitatively 
wrong conclusions about the flow. 

4.2 Dopant Segregation 

Dopant fields were calculated for the velocity field and melt/solid 
Interface shapes discussed above. The results reported here were calculated 
using a 18 x 34 mesh (2553 unknowns) embedded In the 8 x 16 grid used In the 
melt for computing the velocity field. The accuracy of the calculations 
of dopant concentration was assessed by computing the dopant mass balance 
for the entire melt; this balance closed to within five percent for Rayleigh 
numbers up to 1 x 10 . The finite element grids used here were Inadequate 
for dopant calculations at higher values of Ra. 

Dopant fields computed for the segregation coefficient (k » 0.1) and 
Schmidt number (Sc ■ 10) similar to gallium In germanium are shown In Fig, 8. 
The almost parallel Iso-concentration lines for Rayleigh numbers up to 10 
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correspond to the one-dimensional solidification model. The concentration 

field was deformed at higher values of Ra; flow downward along the axis 

swept dopant from the crystal and Induced radial segregation across the 

4 

surface of the crystal. The concentration field at Ra ■ 1 x 1C had the begin- 
nings of the uniform core of melt and steep concentration gradients along 
each boundary which Is consistent with the boundary- layer model for a 
well-mixed melt. A large amount of radial segregation was present even 
at this level of convection. 

The radial variation of dopant composition at the Interface Is shown 
on Figure 9 for the concentration fields In Figure 7. The change In 
concentration across the Interface evolved with Increasing Ra from the 
approximately one percent segregation caused by Interface curvature at Ra * 0 
to almost seventy percent of the nwan value 1/k for the worst case of Ra ■ 1 x 10^. 
The radial variation In dopant decreased with the more Intense flow motion 

4 

that corresponded to Ra * 10 . The maximum In radial segregation with 
increasing convection found here Is consistent with calculations of Nlkltan 
et al. [14] for the growth of gallium-doped germanium In a horizontal boat. 

The percent radial segregation, defined as 

IL £ |c(0,h(0)) - c(A,h(A))| x 100 , (21) 

Is plotted as a function of Ra on Figure 10 for segregation coefficients 
between 0.01 and 0.9. In all cases, Ac reached a maximum for the flow 

3 

corresponding to nearly Ra ■ 10 and decreased for larger values of Rayleigh 
number. The value of the segregation coefficient set the mean concentration 
In the melt at the Interface and affected the level of radial segregation 
In the crystal. The exponentially decreasing concentration profile present 
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at low values of Ra extended further Into the melt for smaller segregation 
coefficients and made the concentration field more sensitive to convection. 
Consequently the level of segregation Ac was higher for lower k. 

Changing the dlffuslvlty of solute also had a marked effect on the 
level of radial segregation, as Is shown on Figure 11 by calculations of 
Ac for a range of Schmidt numbers and all other parameters the same as 
for the gallium-doped germanium system. Increasing Sc decreased the 
length of the concentration gradient at the Interface and caused higher 
Pedet numbers (Pe M = Ra^ 2 Sc/Pr) at lower values of Ra. The concentration 
fields for sample combinations of Sc and P.a are shown In Figure 12. At high 
Schmidt numbers and Rayleigh numbers of 1 x 10^ and above, the Isoconcentra- 
tlon curves develop fingers oriented parallel to the melt/sol Id Interface 
that correspond to extremely rapid variation In composition within a distance 
of the same order of magnitude at the thlcknes* of the diffusion-layer for 
(Ra ■ 0). This variation lead to large (over 100«) radial segregation. 

4.3 Calculation of Effective Segregation Coefficients 

The complicated segregation patterns discussed In the previous section 
create doubts as to the applicability of the simple "diffusion-layer" or 
"boundary-layer" model fcr correlating axial segregation behavior. Our detailed 
results for the velocity and dopant profiles put us In the rare situation 
of being able to check this correlation. To do this, radially averaged 
concentration profiles, defined as 
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c(z) 2 c(r,z) rdr , (22) 

y o 

were calculated and plotted for each set of Ra, k, and Sc; sample profiles 
are shown In Figures 13 and 14 for Sc ■ 10 and Sc ■ 50, respectively. For 
low Rayleigh numbers (0 <, R <_ 1 x 10). these radially averaged profiles are 
essentially the same as the concentration profiles predicted by one-dimensional 
models which account only for the convection caused by crystal growth. A 
region of nearly uniform average concentration c"(z) • c b developed with 

4 

Increasing Ra, as shown on Figure 13; here, the profile for Ra « 1 x 10 can 
be divided Into a region of steep dopant gradient adjacent to the crystal, 
a zone of uniform concentration and a gradient zone caused by the fictitious 
inlet condition at z - h(0) ■ 0. 

The effective segregation coefficient appropriate for a true unsteady 
Bridgman system with A * 0.25 was approximated as 

k eff * c s /e b s 1/e b • (23 > 

where c $ 1$ the dimensionless average concentration across the crystal. Values 
of c b and k fi ^ computed from the finite element results presented 4 n the last 
section are summarized In Table IV. 

If a diffusion- layer of thickness £ 2 6/L exists that separates the 
crystal surface from bulk melt at concentration c b , the effective segregation 
coefficient k 0 ^ can be derived following Flemings [38] as 

* — • . (24) 

* c b k ♦ (l-kJexpf-V^/C) 


or In the dimensionless form used here as 
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k. ff • . (25) 

eTT k + (1-k)erp(-c6 P« Sc/Pr) 

Rearranging eq. (25) yields an explicit relation for the diffusion-layer 
thickness fi: 

s . .SL. ln ° ~ *•"! I • ( 26 ) 

PeSco I k eff (1-k) ) 

Values of 6 computed from eq. (25) with k g ^ calculated from the finite 
element results are compared In Table IV to values 6^ measured from the 
profiles of F(z). In each case d and d^ are of the same order of magnitude 
but may differ by as much as a factor of two. This result Is not surprising; 
within the distance 6 f from the crystal In the finite element calculations, 
the velocity field Is a combination of the growth velocity and a contribution 
from natural convection. These two components are of the same order of 
magnitude and It Is the latter component which causes the discrepancy between 
6 and 6^ and which leads to radial segregation. 

4.4 Variation of Thermal Conductivities: K 5 kj/k^ 

The ratio of thermal conductivities K s kj/k^ Is se"dom unity for real 
materials, as has been assumed In the calculations presented above. Metallic 
semiconductors usually have higher conductivities In the melt (K < 1), while 
oxide systems usually have greater conductivities In the solid phase (K > 1). 
We have examined the effect of varying the conductivity ratio by changing K 
to 0.5 and 2.0. Temperature fields are shown In Figure 15 for these two cases 

4 

at Ra ■ 1 x 10 and are essentially the same as for lower values of Rayleigh 
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number. Changing K from unity destroyed the axial synmetry In the temperature 
Meld and displaced the Interface from the center of the adiabatic zone. For 
K ■ 0.5, the melt penetrated deeper Into the ampoule and the Interface was 
convex toward the crystal. The Interface was concave for K ■ 2.0 and located 
above the center of the adiabatic zone. In each case, the deflections of the 
Interface were nearly ten percent of the mean location compared to the four 
percent deflections calculated when the conductivity ratio was unity. 

Flow fields are displayed In Figure 16 for K ■ 0.5 and 2.0 and Rayleigh 

4 

numbers between zero and 1 x 10 . The changes In conductivities have little 
effect on the flow pattern besides the small deflection In the flow near the 
solidification Interface needed because of the higher curvature of the crystal 
surface. The concentration fields shown on Figure 17 are much more distorted 
by the Interface curvature; at low convection levels (0 ^ Ra 1 x 10 ) significant 
radial segregation occurred, as marked by the Intersection of Iso-concentration 
curves with the Interface. The higher Interfaclal dopant concentrations appeared 
where the Interface had Its largest positive curvature as measured from the 
melt; this point shifted from the centerline for K ■. 0.5 to the wall of the 
ampoule for K ■ 2.0. The amount of radial segregation Ac caused solely 
by Interface curvature was 23t for K ■ 0.5 and 18S for K ■ 2.0. 

Fluid motion caused by natural convection swept dopant from the center 
of the Interface towards the ampoule wall. For K ■ 0.5, weak convection 
decreased the radial segregation by stripping away the dopant peak caused 
by Interfaclal curvature. Cellular convection In the melt Increased the radial 
segregation for the case with K ■ 2.0 by adding to the layer of dopant 
already at the edge of the crystal. These results are suwnarlzed on 
Figure 18 as a plot of Ac as a function of Ra. 
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4.5 Variation of Malt and Adiabatic Zone Lengths 

The radial temperature gradients that drive natural convection In the 
vertically stabilized Bridgman configuration are not physically modelled 

3 

by the Rayleigh number Ra s $gl (T h - T c )/a & v scaled with the length of the 
ampoule. In fact, the structure of the flow should be Independent of ampoule 
length as long as the flow cell does not Interact with the top (z ■ 0) of 
the ampoule. The size of the radial temperature gradients will be propor- 
tlonal to that of the axial gradient In the gradient zone (T h - 1 C )/L 
and so a more physically correct measure of the buoyancy driving force Is 



Calculations were performed for melt lengths of 2R, 3R and 4R, with the crystal 
length and all other dimensionless groups held at the values of the base case 
studied above; the overall ampoule lengths L were 4R, 5R and 6R, respectively. 
The Intensity of the convection as measured by Ra f was held constant for each 
aspect ratio of the melt. Temperature , velocity and dopant fields are shown 
In Figure 19 for Ra r ■ 1.5625, which corresponds to Ra ■ 1 x 10 3 for the base 
case A ■ 0.25. It Is clearly seen that these three fields are essential Ty 
Identical for all three aspect ratios. The aspect ratio A ■ 0.25 is adequate 
for modelling much longer ampoules. 

Tha calculations presented In Sections 4.1 - 4.3 show that radial 

segregation results from the Interaction of the cellular flow driven by the 

radial temperature gradients at the edge of the adiabatic zone with the 

concentration gradient adjacent to the melt/sol Id solid Interface. 

* 

Hence, the length L g z L g /L of the adiabatic zone should be an Important 
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variable In determining the extent of this segregation. Calculations were 
performed for an ampoule with aspect ratio A ■ 1/8 and adiabatic zones with 
dimensions L g of 1/8 (same as In Sections 4.1 - 4.3), 1/4, 3/8, and 1/2. 

Samples Isotherms, streamlines and concentration fields are shown In 

3 

Figure 20 for Ra ■ 8.0 x 10 and the thermophyslcal properties corresponding 
to the gallium-doped germanium system listed In Table II. As the length 
of the gradient region was Increased the center of the cellular flow moved 
with the edge of the adiabatic region and was located farther from the melt/ 
solid Interface. For the longest gradient zone l g ■ 1/2 (see Fig. 20d) the 
bulk of the convection occurred beyond the concentration gradient emanating 
from the melt/sol Id Interface and the concentration field there was only 
slightly disturbed by the flow. Very little radial segregation Is expected 
for these conditions. The percent radial searegatlon Ac Is plotted on Fig. 20 
as a function of Rayleigh number fur the four values of L g listed above. 

Fixing Ra and Increasing L g corresponds experimentally to holding 
constant the overall temperatu.e difference (T h - T c ) and changing only 
the length of the adiabatic zone. This procedure decreases the temperature 
gradient In the melt and weakens the flows, as shown on Fig. 20 by the 
increased spacing of the Isotherms and by the decrease of the circulation 
*max Increasing L g . Calculations with constant Intensity of the 
convection require fixing the radially-scaled Rayleigh number Ra r (see 
eq. 27). Curves of constant Ra r are shown on Fig. 21 and demonstrate 
tht' decrease In radial segregation with Increasing L g . At moderate convection 
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level s the adiabatic zone can bt nada long tnouqh that the cellular flow 
does not Interact with the dopant diffusion-layer adjacent to the crystal. 

The minimum length for this coupling will he strongly dependent on the value 
of the dlffuslvlty (Schmidt number), which sets the length of the concentra- 
tion gradient; see Fig. 10. 
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5. VERTICALLY DESTABILIZED CONFIGURATION 

When the melt Is positioned below the crystal the vertical temperature 
gradient leads to unstable density stratification which enhances convection. 
Calculations were performed for the destabilized geometry over a range of 
Rayleigh number Ra with the the rmophys leal properties listed In Table II. 

The size of the ampoule was taken to be the same as the base case for 
the stabilized system, l.e. L g ■ 1.0 and A * 1/4. For Rayleigh numbers up 
to Ra - 1 x 10 5 , the temperature fields for the destabilized configuration 
were indistinguishable from the one for pure conduction (Ra ■ 0; see Fig. 3a) 
and.: so are not shown here. Streamlines computed for four representative 
values of Ra are displayed as Fig. 22. 

The direction of the How In the destabilized system was opposite that 
for the stabilized configuration *nd had qualitatively different effects on 
melt/solid Interface shape and solute segregation. Melt moved upward along 
the sidewall of the ampoule, turned Inward at the melt/solld Interface, 
and fell along the centerline of the melt. As In the stabilized geometry, 
three regimes of the flow were located with changing Ra. Cellular fluid 
motion developed near the centerline for low values of Rayleigh number 
and spread toward the ampoule wall with Increasing Ra. The center of the 
cell was positioned near the edge of the adiabatic zone. Indicating that 
the radial temperature gradients at this edge were Influencing the flow. 

At Ra ■ 1 x 10 5 (see Fig. 22d) the destab'llzlng axial temperature gradient 
was controlling the flow and the cell had migrated downward In the ampoule. 
Calculations In a similar flow geometry and low Prandtl number [31] have 
indicated that the destabilized flow will separate from the sidewall at 
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higher Rayleigh numbers. 

Melt/solld Interface shapes are shown In Fig. 23 for the flow fields 
given In Figure 22. Again, conduction dominated the heat flux In the melt 
near the surface of the crystal and the shape of the solidification Isotherm 
was only slightly altered by convection. The shape shown for Ra between 0 
and 1 x 10^ was the same as shown for the stabilized configuration and low 
Rayleigh numbers In Fig. 6. The hot melt rising along the ampoule wall 

3 

pushed the Interface higher Into the adiabatic zone at Ra * 1 x 10 and 
above and resulted In Interfaces with convection being "flatter" than the 
shape corresondlng to Ra • 0. The maximum height of the Interface for Ra * 

1 x 1(T occurred off the centerline of the ampoule. 

More dramatic differences are found by comparing radial dopant profiles 
for the two geometries. Dopant fields and Interface distributions are shown 
In Figs. 2d and 25 for the parameters listed In Table II. The flow upward 
to the Interface swept dopant from the outer edges of the crystal to the 
center of the ampoule. This led to Interfaclal concentrations that were 
highest at the center of the crystal, see Fig. 24. The amount of radial 

4 

segregation across the crystal decreased for Ra ■ 1 x 10 , just as In the 
case of the stabilized configuration. The percent radial segregation expected 
In the stabilized and destabilized geometries are compared directly on Fig. 26. 
The maximum In radial segregation Is twenty percent less for the destabilizing 
configuration than for the stabilized system. 
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6. DISCUSSION 

The coupled analysis of natural convection, irelt/solid interface shape 
and dopant segregation In small-scale crystal growth systems is complicated, 
but Is no longer impossible, because of the advent of new finite-element 
methods. The detailed pictures of the field variables and solidification 
front for the vertical Rridgman system studied here give new insights into 
the role of convection In setting the effective segregation between the bulk 
melt and crystal and the radial segregation in the crystal. 

For moderate levels of convection# the stagnant-film or diffusion-layer 
model is a gross over-simplification of these Interactions. Diffusion-layer 
thicknesses 6 determined from eo. (26) and experimental data are, at best, 
an empirical fit to the effective segregation coefficient for the crystal/melt 
system. The comparison In Section 4.3 demonstrates that, although the radially 
averaged diffusion-layer thickness determined from the finite element simulations 
and eg. (26) are of the same order-of-iragnitude, the actual concentration gradient 
adjacent to the crystal is far from radially uniform. As much as sixty percent 
radial segregation can exist. Only detailed calculations of the exact Interaction 
of fluid flow and the dopant profile adjacent to the Interface can be used 
to estimate the level of radial segregation in the crystal. 

Our study of a prototype Bridgman system gives qualitative understanding 
of fluid flow and dopant segregation In actual growth systems and will serve 
as the starting point for more refined calculations aimed directly at comparison 
with experiments In a well -characterized small-scale system [36]. In these 
calculations, the ratio of thermal conductivities K between melt and crystal 
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and the length of the gradient zone L q have emerged as important design 

parameters for determining Interface shape and radial segregation. The 

value of K sets the shape of the Interface and determines the Qualitative 

behavior of the radial segregation with Increased convection. The length 

L„. sets the degree of Interaction between the cellular flow and the 
S 

concei. ..ration gradient adjacent to the Interface. For long enough gradient 
zones and moderate convection levels, the flow causes little radial segrega- 
tion; see Fig. 21. 

Accurate comparisons with experimentally determined interface shapes 
and radial and axial dopant profiles will require the Inclusion of the 
temperature variations of key thermcphyslcal properties, such as the melt 
and solid conductivities and the density of the melt. Most of these variations 
are unknown. The thermal Interactions between the ampoule and the surroundings, 
e.g. heat capacity and conductivity of the ampoule, will also have to be 
more accurately accounted for to correctly model radial temperature gradients 
In the system. 

The slow translation of the ampoule has little effect on the temperature 
and velocity fields, when the ampoule Is sufficiently long; however, the average 
concentration of dopant In the melt will Increase slowly with time. For con- 
vection levels Intense enough that a plateau In c“(z) forms In the melt (see 
Fig. 13), this transient can be studied by an unsteady dopant calculation 
with fixed melt length and steady fluid flow, but with a time varying Inlet 
concentration determined by an overall mass balance. For very low convection 
levels and high dlffuslvltles (low Sc) where the concentration gradient extends 
to the end of the ampoule, fully transient calculations taking Into account 
the length of the melt will be necessary. Efficient numerical methods presen +, y 
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exist only for one-dimensional solidification problems [39] and are being 
extended to two-dimensional simulations. 

The most critical issue not addressed in this paper Is the transition 
from the steady-state flow patterns predicted here to time-periodic flows 
observed In some experiments [5, 29]. Mathematically, these transitions 
occur as Hopf bifurcations from the steady flows and have been predicted 
theoretically by computer-aided bifurcation analysis combined with finite- 
element calculations [40] for low Prandtl number melts in the horizontal 
boat geometry studied experimentally by Hurle et al. [41], Simulation 
of crystal growth In the time-periodic regime requires time-dependent 
solution of the Bouslnesq equations as performed by some researchers 
[14, 42] without the Interaction with the melt/solid Interface. Fully 
transient calculations are no panacea. Besides belna extremely expensive, 
their accuracy Is difficult to assess. The loss-of-exlstence of steady 
solutions shown In Section 4.1 for crude finite element approximations 
would be manifested as time-fluctuating flows In transient calculations. 

As is the case here, refined transient calculations may lead to steady 
velocity and temperature fields. 
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FIGURE CAPTIONS 


Dimensionless Parameters and Representative Values 
for the Gallium-Doped Germanium System. 

Material Properties Characteristic of Gallium-Doped 
Germanium. 

Accuracy of Finite Element Calculations of Flow and 
Temperature Fields. 

Comparison of Effective Segregation Coefficients 
Predicted by Finite Element Analysis with Diffusion- 
Layer Model, Eq. (26). 

Geometry of prototype vertical Bridgman growth system. 

Sample finite element mesh for computing velocity and 
temperature fields and melt/solld Interface shape. 

Sample flow fields for vertically stabilized Brldernmn 
system with Rayleigh numbers between 0 and 2.5 x 10 b . 
Streamlines are spaced at equal Intervals between the 
maximum (or minimum) value for the cell and zero. 

Maximum streamfunctlon iw,ax as a function of Ra for 
vertically stabllfzed Bridgman system. 

Sample temperature fields for the vertically stabilized g 
Bridgman system with Raylelqh numbers between 0 and 2.5 x 10° 

Melt/solld Interface shapes for the flow and temperature 
fields shown In Figs. 3 and 5. 

Change of structure with Rayleigh number of flow families for 
three finite element meshes, as represented by the 
Nusselt number Nu^. 

Dopant concentration fields for flow fields shown In 
Fig. 3; k ■ 0.1 and Sc • 10. 

Variation of Interfaclal dopant concentration with 
Rayleigh number, k ■ 0.1 and Sc ■ 10. 

Percent radial segregation as a function of Rayleigh 
number and segregation coefficient for the vertically 
stabilized Bridgman system; Sc ■ 10. 
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Figure 21. 

Figure 22. 

Figure 23. 
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Percent radial segregation as a function of 
Rayleigh number and Schmidt number for the 
vertically stabilized Bridgman system; k ■ 0.1. 

Dopant concentration fields for three Schmidt 
numbers as a function of Rayleigh number; k * 0.1. 

Profiles of radially averaged dopant concentration 
for three values of Ra; k ■ 0.1 and Sc * 10. 

Profiles of radially averaged dopant concentration 
for three values of Ra; k ■ 0.1 and Sc ■ 50. 

Sample temperature fields for vertically stabilized 
Bridgman systems with Ra ■ 1 x 10 4 and thermal conduc- 
tivity ratios of (a) K ■ 2.0 and (b) K ■ 0.5. 

All other parameters are the same as listed In Table I. 

Sample flow fields for vertically stabilized Bridgman 
systems with thermal conductivity ratios of (a-d) 

K ■ 0.5 and (e-f) K - 2.0. 

Dopant concentration fields for flow fields shown In 
Figure 16; k » 0.1 and Sc ■ 10. 

Percent radial segregation as a function of Rayleigh 
number and thermal conductivity ratio; k ■ 0.1 and 
Sc ■ 10. 

Effect of melt length on temperature, velocity and 

dopant fields for Ra r ■ 1.5625 * k • 0.1 and Sc ■ 10. 

Effect of the length of the gradient zone l. on 
temperature, velocity and dopant fields for q 
Ra ■ 1 x 104, k * 0.1 and Sc ■ 10. 

Effect of the length of the gradient zone L on 
the percent radial segregation Ac; k ■ 0.'! and 
Sc ■ 10. 

Sample flow fields for vertically destabilized 
Bridgman system with the thermophysical properties 
given In Tables I and II. 

Melt/sol Id Interface shapes for the flow fields 
shown In Fig. 22. 

Dopant concentration fields for flow fields shown 
In Fig. 22; k ■ 0.1 and Sc ■ 10. 


Figure 25 
Figure 21 
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Variation of Interfaclal dopant concentration 
with Rayleigh number; k ■ 0.1 and Sc * 10. 

Comparison of percent radial segregation for 
stabilized and destabilized Bridgman systems; 
k ■ 0.1 and Sc ■ 10. 
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Table I. Dimensionless Parameters and Representative 
Values for Gallium-Doped Germanium System 


GROUP 

DEFINITION 

VALUE 

Rayleigh Number 

Ra s B<](T h - T e )L 3 /« 4 v i 

0 - 10 7 

Prandtl Number 

Pr s v/a £ 

0.01 

Peclet Number 

P. a V,L/a t 

0.01 

Conductivity Ratio 

K i k,/k t 

1.0 

Stefan Number 

S 5 4H f /P t C pi <T h . T e ) 

1.0 

Schmidt Number 

Sc = v/O 

10 

Thermal Dlffuslvlty 
Ratio 

Y 5 0 $ /0 t 

1.0 

Density Ratio 

0 s 0j /o t 

1.0 
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Table II. Material Properties Characteristic of Gal Hum-Doped Germanium 


Property 

Thermal conductivity of melt (k^) 
Thermal conductivity of solid (k s ) 
Heat capacity of melt (c p ^) 

Heat capacity of solid (c D ) 

H s 

Density of melt (p^) 

Oenslty of solid (p s ) 

Melting temperature (T^) 

Kinematic Viscosity of Melt (v) 

Heat of fusion (AH { } 

Thermal Expansion Coefficient (6) 

Dlffuslvlty of Ga In Ge (fl) 

Segregation Coefficient of 
Ga In Ge (k) 


Value 

0.17 W/K cm 
0.17 W/K cm 
0.17 J/gm K 
0.17 J/gm K 
5.6 o/cm 3 
5.6 g/cm 3 
956*C 

2 

0.13 cm /sec 
506 J/gm 

0.25 x 10‘ 3 (Kf 1 
1.3 x 1 O’ 4 cm 2 /sec 


0.1 
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MAXIMUM OF STREAMFUNCTION.Y. 
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Figure Ar 
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AXIAL COORDINATE 
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INTERFACE NUSSELT NUMBER Nu^ (10 
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Figure 7 
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CONCENTRATION PROFILE AT INTERFACE c(r,h<r)) 
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DIMENSIONLESS RADIUS, r 
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RADIALLY AVERAGED DOPANT CONCENTRATION c(z) 





RADIALLY AVERAGED DOPANT CONCENTRATION cU) 
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DIMENSIONLESS DISTANCE FROM INTERFACE.z-h (r) 


Figure 14 






OF POOR 


(o) K* 2.0 


(b)K«0.5 



Pi C^ uvc. 15 

CUaMa au^\ £>\focD*i 

















PERCENTAGE RADIAL SEGREGATION, AC * 100% 
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PERCENTAGE RADIAL SEGREGATION AC x 100% 
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(a) Ro*0 


(b)Ra*lilO s (c)Ro»UI0 4 



(d)Ro>tilO s 
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RADIAL COORDINATE, r 
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PERCENTAGE RADIAL SEGREGATION , AC 
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